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Thermal evolution of a pulsating neutron star 
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ABSTRACT 

We have derived a set of equations to describe the thermal evolution of a neutron 
star which undergoes small-amplitude radial pulsations. We have taken into account, 
in the frame of the General Theory of Relativity, the pulsation damping due to the 
bulk and shear viscosity and the accompanying heating of the star. The neutrino 
emission of a pulsating non-superfluid star and its heating due to the bulk viscosity 
are calculated assuming that both processes are determined by the non-equilibrium 
modified Urea process. Analytical and numerical solutions to the set of equations of the 
stellar evolution are obtained for linear and strongly non-linear deviations from beta- 
equilibrium. It is shown that a pulsating star may be heated to very high temperatures, 
while the pulsations damp very slowly with time as long as the damping is determined 
by the bulk viscosity (a power law damping during 100-1000 years). The contribution 
of the shear viscosity to the damping becomes important in a rather cool star with a 
low pulsation energy. 

Key words: stars: neutron - evolution - oscillations. 



1 INTRODUCTION 

Dissipation processes play an important role in the neu- 
tron star physics; for instance, they determine the damping 
of stellar pulsations (see, e.g., Cutler, Lindblom & Splinter 
1990). Pulsations may be excited during the star formation 
or during its evolution under the action of external pertur- 
bations or internal instabilities. The instabilities arising in a 
rotating star from the emission of gravitational waves may 
be suppressed by dissipation processes. This affects the max- 
imum rotation frequency of neutron stars and creates prob- 
lems in the detection of gravitational waves (see, e.g., Zdunik 
1996; Lindblom 2001; Andersson & Kokkotas 2001; Arras et 
al. 2003). 

The joint thermal and pulsational evolution of neutron 
stars was studied long ago (e.g., Finzi & Wolf 1968 and ref- 
erences therein). Naturally, it was done with a simplified 
physics input and under restricted conditions (Section^. 
However later, while estimating the characteristic times of 
pulsational damping, one usually ignored the temporal evo- 
lution of the stellar temperature (see, e.g., Cutler & Lind- 
blom 1987; Cutler et al. 1990), which led to an exponential 
damping. This is not always justified because the parameters 
defining the damping rate, e.g., the bulk and shear viscosity 
coefficients, are themselves temperature-dependent. 



Clearly, the temperature variation can be neglected if 
the characteristic damping time r <C t coo i and -Epuis <C Eth, 
where i coo i is the characteristic time of neutron star cooling, 
while -Epuis and E t h are the pulsational and thermal ener- 
gies, respectively. We will show that these conditions are 
violated in a wide range of initial temperatures and pulsa- 
tion amplitudes. 

This paper presents a self-consistent calculation of the 
dissipation of radial pulsations with account for the ther- 
mal evolution of a non-superfluid neutron star whose core 
consists of neutrons (n), protons (p) and electrons (e). It 
extends the consideration by Finzi & Wolf (1968) (see Sec- 
tion^Jfor details). We consider two dissipation mechanisms: 
one is via the non-linear (in the pulsation amplitude) bulk 
viscosity in the stellar core and the other is due to the shear 
viscosity. We neglect other possible dissipation mechanisms, 
particularly, the damping of pulsations induced by the star 
magnetic field (as discussed in detail by McDermott et al. 
1984 and by McDermott, van Horn & Hansen 1988). The 
magnetic field is assumed to be low. 



2 EIGENFUNCTIONS AND 
EIGENFREQUENCIES OF 
NON-DISSIPATIVE RADIAL PULSATIONS 



* E-mail: gusakov@astro.ioffe.ru (MEG); yak@astro.ioffe.ru 
(DGY); ognedin@astronomy.ohio-state.edu (OYG) 



Here we discuss briefly radial pulsations of a neutron star, ig- 
noring energy dissipation. This problem was first considered 
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by Chandrasekhar (1964), and we will refer to his results. 
The metric for a spherically symmetric star, which experi- 
ences radial pulsations, can be written as 

ds 2 = -e v dt 2 + r 2 dQ. 2 + e A dr 2 , (1) 

where r and t are the radial and time coordinates and dQ. 
is a solid angle element in a spherical frame with the origin 
at the stellar center. Here and below, we use the system of 
units, in which light velocity c = 1. The functions v and 
A depend only on r and t and can be written as: v(r,t) = 
vo(r)+8v(r, t); X(r,t) = Xo(r)+8X(r,t). Herefo(r) and Ao(V) 
are the metric functions for an unperturbed (equilibrium) 
star, and 5v(r,t) and SX(r, i) are the metric perturbations 
due to the radial pulsations (described by Eqs. (36) and (40) 
of Chandrasekhar 1964). 

The radial pulsations can be found by solving the 
Sturm-Liouville problem (Eq. (59) of Chandrasekhar 1964). 
A solution gives eigenfrequencies of pulsations uik and eigen- 
functions fjfe(r), where ffc(r) is the Lagrangian displace- 
ment of a fluid element with a radial coordinate r. By 
neglecting the dissipation and the non-linear interaction 
between the modes (the pulsation amplitude is taken to 
be small, |£fe(r)| <C r), we can write the general solution 
for a fc-th mode as £,(r,t) = £fc(r) cosuikt- The boundary 
conditions for the Sturm-Liouville problem have the form: 
P(r = R + £(R,t)) = 0, £(0,t) = 0, where P(r,t) is the 
pressure and R is the unperturbed stellar radius. 

We employ the equation of state of Negele & Vautherin 
(1973) in the stellar crust and the equation of state of Heisel- 
berg & Hjorth- Jensen (1999) in the stellar core. The lat- 
ter equation of state is a convenient analytical approxima- 
tion of the equation of state of Akmal & Pandharipande 
(1997). For this equation of state, the most massive sta- 
ble neutron star has the central density p c = 2.76 x 10 15 
g cm -3 , the circumferential radius R — 10.3 km, and the 
mass M = M max = 1.92 Mq. The powerful direct Urea pro- 
cess of neutrino emission is open in the core of a star of mass 
M > 1.83 M Q . 

An important parameter which enters the equation of 
radial pulsations is the adiabatic index 7. Since the fre- 
quency of stellar pulsations is u>k 3> 1/turca, where t\j Ica is 
the characteristic beta-equilibration time (see, e.g., Haensel, 
Levenfish & Yakovlev 2001; Yakovlev et al. 2001), the adi- 
abatic index must be determined assuming the "frozen" 
nuclear composition (see, e.g., Bardeen, Thorne & Meltzer 
1966): 

din P(n h ,x c ) 
7= dlnn b ' (2) 

where rib is the baryon number density, x c — n c /rih, and n c 
is the electron number density. 

The relative radial displacement of matter elements in 
a pulsating star (in the absence of dissipation effects) will 
be described by a small parameter e: 

e = lim f fc (r)/r. (3) 

r— >0 

Thus, e determines the normalization of the function £fc(r). 

Figure 1 shows the dependence of £k(r)/r (artificially 
normalized such that |e| = 1) on the distance to the stellar 
center r for the first three modes with the frequencies u)a = 
1.705 x 10 4 s" 1 (solid line), Wi = 4.121 x 10 4 s -1 (long 
dashed line), and U2 = 5.950 x 10 4 s _1 (short dashed line), 



respectively. By way of illustration, we consider a model of 
a star of mass M = 1.4 M© (R = 12.17 km, p c = 9.26 x 10 14 
g cm -3 ). As expected, the fundamental mode is close to the 
homological solution £o(r) = r. Introducing a normalization 
constant, we get fo(r) = sr. Therefore, for the fundamental 
mode, e determines the amplitude of relative displacements 
of the pulsating stellar surface. 

We will further need the pulsation energy, which can 
be calculated if we formulate, for example, the variational 
principle for the characteristic eigenvalue problem in ques- 
tion. For the fc-th radial mode, we have (see, e.g., Meltzer & 
Thorne 1966) 

£puis = \ J (P + P) [e {X °~ va)/2 ^ &] 2 e"°/ 2 dV , (4) 

where p is the mass density and dV = 47rr 2 e Ao//2 dr is the 
volume element measured in a comoving frame. 

The account of the energy dissipation in the fc-th mode 
leads to a relatively slow damping of pulsations. In particu- 
lar, we take for the Lagrangian displacement 

£(r,t) = C fc (*)&(r) coswfct. (5) 

where Ck(t) is a slowly decreasing function of time (the char- 
acteristic dissipation time r 3> 1/uJk), which will be further 
termed the pulsation amplitude. The dissipation is assumed 
to be "switched on" at the moment of time t = 0, at which 
the initial amplitude is 

C fc (0) = 1. (6) 

From Eq. Q the pulsation energy in the fc-th mode with 
dissipation is 

-Epuis(t) = -Epuiso Cfc(t). (7) 

Using Eqs. Q and Q, we can estimate the pulsa- 
tion energy for the fundamental mode, E pu i s (t) ~ 2 x 
10 53 L0 2 e 2 Cl(t) erg. The thermal energy of the star is E th ~ 
(47r/3) R 3 ct T ~ 10 48 T 9 2 erg, where ct oc T is the spe- 
cific (per unit volume) heat capacity of the stellar matter 
(see, e.g., Yakovlev, Levenfish & Shibanov 1999) and T9 
is the internal temperature of the star in units of 10 9 K, 
(x)4 = LUk/{W 4 s -1 ). These estimates show that there is a 
wide range of values of the parameters e, C'k(t), and T, at 
which -Ep U i s > -Eth- In such a case, one should account, at 
least, for the stellar temperature evolution during the damp- 
ing of pulsations. 



3 NON-EQUILIBRIUM MODIFIED URCA 
PROCESS 

The condition for beta-equilibrium in the stellar core has 
the form: Sp(r, t) — p n — p p — p c = 0, where pi is the chemi- 
cal potential of particle species i = n, p, e. Stellar pulsations 
lead to deviations from beta-equilibrium, &p(r, i) 7^ 0, and 
hence, to the dissipation of the pulsation energy. The dissi- 
pation rate is determined by the processes tending to return 
the system to equilibrium. We suggest that the direct Urea 
process in the neutron star core is forbidden. Then the main 
process which determines the pulsation energy dissipation is 
the modified Urea process. In this section, we will discuss the 
non-equilibrium modified Urea process and obtain the rela- 
tionship between the Lagrangian displacement £,(r,t) and 
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Figure 1. The parameter £k/ r normalized such that [e| = 1, for the fundamental, first, and second modes of radial stellar pulsations 
(solid, long-dashed, and short-dashed lines, respectively) versus the dimensionless radial coordinate r/R. 



the parameter S/j,(r, i) that characterizes the local deviation 
from beta-equilibrium. 

The non-equilibrium modified Urea process has been 
discussed since the end of the 1960s (see the classical pa- 
per by Finzi & Wolf 1968 and references therein). How- 
ever those old studies were accurate only qualitatively (see, 
e.g., Haensel 1992 and Section Later the problem has 
been reconsidered by several authors (see, e.g., Haensel 1992; 
Reisenegger 1995; Haensel et al. 2001). Below, references will 
primarily be made to the review of Yakovlev et al. (2001) 
since we employ similar notations. The modified Urea pro- 
cess has the neutron and the proton branch, each including 
a direct and an inverse reaction: 

n + N^p + N + e + ^o, p + N + e^n + N + ^ c . (8) 

Here N = n or p for the neutron or proton branch, respec- 
tively. In beta-equilibrium, the neutrino emissivities in these 
two channels are given by 

qM „ g.ixW^W^P) 
\m n j \m p J 

( \ 1//3 

x f-£J TfctuP* erg cm" 3 s" 1 , (9) 

Q$ - «£? r4V ^ + >-^ )2 e, (io) 



where no = 0.16 fm -3 is the nucleon number density in 
atomic nuclei; n p is the proton number density; m n and m p 
axe the masses of free neutrons and protons; m„ and m P are 
the effective masses of neutrons and protons in dense matter; 
p Fe ,p Fp , and p Fn are, respectively, the Fermi momenta of 
electrons, protons, and neutrons; and a n ,f3 n ~ 1 are the 
correction factors (for details, see Yakovlev et al. 2001). In 
Eq. 11011 the function 0=1 if the proton branch is allowed 
by momentum conservation (p^ < Sppp +p Pc ), and O = 
otherwise. 

In beta-equilibrium, the direct and inverse reaction 
rates in Eq. @ coincide, i.e., the matter composition does 
not change with time. The reactions involve only parti- 
cles in the vicinity of their Fermi surfaces, with the energy 
\ e i ~ Mi| ksT, where i — n,p,e, and fee is the Boltzmann 
constant. Therefore, the neutrino emissivity depends sensi- 
tively on the temperature, and the process cannot occur at 
T — 0. A drastically different situation arises in the presence 
of deviations from beta-equilibrium (<5u 7^ 0). The rates of 
the direct and inverse reactions become different, the system 
tends to equilibrium, and the matter composition changes; 
the process remains open even at T — 0. 

Let r and V be the numbers of direct and inverse re- 
actions of the modified Urea process per unit volume per 
unit time. The analytical expressions for Ar = T — F and 
for the neutrino emissivity Qnoneq = + Q^ p ' of the non- 
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equilibrium modified Urea process were derived by Reiseneg- 
ger (1995): 



Ar 



Qn 



11513 k B T 
Q cq F(y), 



yH(y), 



(11) 
(12) 



where Q eq = 
are given by 



+ Qcq', and the functions H(y) and F(y) 



Icq 



H(y) 
F{y) 



189tt 2 y 21tt 4 t/ 4 
1 + — T^r- + ■ 



367 



367 



+ 



o 6 6 

1835 ' 



, 22020tt 2 y 1 5670tt 4 y 4 
1 + — „ + 



11513 

,.6 



11513 



+ 



4207r b y° 9tt s y a 
11513 + 11513 ' 



(13) 



(14) 



Here y = 8[i/(-K 2 k B T)\ the factor it 2 in the denominator 
is introduced to emphasize that the real variation scale of 
the functions H(y) and F(y) is Sfi/(10kBT) (but not just 
Sfi/ksT). It follows from Eqs. HI IB and 1121 that there are 
two pulsation regimes. The regime with <5/j <C k^T (y <C 1) 
will be referred to as subthermal and that with S/j, 3> k B T 
(y S> 1) as suprathermal. Relative displacements of fluid 
elements in both regimes are taken to be small (e <C 1). 
From these equations one can see that the values of <5/iAF 
and Qnoncq in the suprathermal regime are independent of 
temperature. 

Let us now find the relationship between the Lagrangian 
displacement £(r, t) and the chemical potential difference 
Sfi(r,t). The quantity Sfi can be treated as a function of 
three thermodynamic variables, say, n b ,n c and T: S/j, = 
Sfi(n b ,n c ,T). During pulsations, these variables will de- 
viate from their equilibrium values ribo, n^o and To by 
An b (r,t), An c (r,t) and AT(r, t). Taking the deviations to 
be small, i.e. obeying the inequality e <C 1, one can write: 



5fi(r, t) 



d5fi(n b0 ,n e0 ,T ) 



+ 



dn b0 
dSfi(n bo ,n e0 ,To) 



dn c 



An e (r,t) 



An b (r,t) 

dSfi(n b0 ,n e0 ,T ) 



dT 



AT(r,t). 
(15) 



The last term in Eq. 11511 can be neglected because 
dSn(n ho ,neo,T )/dTo oc T and AT(r,t) ~ An b (r, t) T /n ba 
(see, e.g., Reisenegger 1995). Accordingly, for a strongly de- 
generate matter (/x, 3> ksTo, i = n,p, e), this term is much 
smaller than the first two terms. The temperature T will fur- 
ther denote an "average" temperature To and its oscillations 
around the equilibrium value will be neglected. 

The form of the functions n b (r, t) and n c (r,t) can be 
found from the continuity equations for baryons and elec- 
trons: 



are quadratic and higher order in f(r, t), one obtains: 

= 0, (18) 
= Are" 072 . (19) 



dn^ + e^ d_ ( bor 2 e -„ /2 dgr,t) 
dt r 2 dr \ dt 



dn e e"° /2 d_ 
dt r 2 dr 



ncorV , o/2 ^) 

at 



These expressions have been derived using Eq. (36) of Chan- 
drasekhar (1964) for the correction SX(r, i) to the metric (see 
Section 2): 



8\[r,t) = -£{r : t) — (A + f ). 
Equation 1181 is easily integrated and yields 



(20) 



An b (r,t) 



n b (r,t) — n b0 
e"°/ 2 d 



-Qr 1 Ubor e 



(21) 



The solution to Eq. 1191 can be written as 
An e (r,t) = n c (r,t) - n c0 = An c0 (r,t) + An cl (r,t), (22) 
An c0 (r,t) = (n c0 r 2 e-^ 2 £(r,t)) , (23) 



where An e o(r, f) describes variations of the electron num- 
ber density ignoring beta-processes, while the function 
Ariel (r,t) describes variations determined by these pro- 
cesses. The latter function satisfies the equation 



OA 



iid 



M 



Are 



(24) 



Generally, the source AF is a complicated function of the 
electron number density n c (r, i). We are, however, interested 
in the high frequency limit, where cut 2> 1/furca (see Section 
2). In that case, the source in the right-hand side of Eq. 
1191 is smaller than other terms. This means that changes 
in the electron number density due to beta-transformations 
are relatively small in a pulsating star (see, e.g., Haensel 
et al. 2001). Therefore, the small parameter An c i can be 
omitted in Eq. 1)15^ . 

By substituting the expressions for An b and An e from 
Eqs. i'2 l\i and l|23(l into Eq. I|150 . we find the relationship 
between Sfi(r,t) and £(r, t): 



= J^Ko,^o) nM a / 2e _, 0/2 

oribo t dr \ 

(25) 

Note that the partial derivative with respect to n b o is taken 
at constant x a o = n c o/n b o- Using Eq. (12511 . we can express 
the parameter y = 5/i/(7r 2 fcBT), as well as AF and Qnoneq 
(see Eqs. 1111 and 1121 1. through the Lagrangian displace- 
ment £(r, i). The relationship between Sfi(r, t) and f (r, i) for 
non-radial pulsations can be derived in a similar way. 



(n b u a ): 
(n c u a ); 



= 0, 

= Ar. 



(16) 
(17) 



Here u a = da; a /ds is the velocity four- vector of the pulsating 
matter. Note that the source AF in the continuity equation 
for electrons is responsible for beta-relaxation processes. 

Writing explicitly the covariant derivatives in Eqs. I|16|l 
and 1171 in the metric Q and neglecting all terms which 



4 THE EQUATIONS OF STELLAR THERMAL 
EVOLUTION AND PULSATION DAMPING 
OUT OF BETA-EQUILIBRIUM 

The thermal balance equation for a pulsating neutron star 
will be derived taking into account three dissipation mech- 
anisms: the shear viscosity in the core, the non-equilibrium 
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beta-processes in the core, and heat conduction. The equa- 
tions of relativistic fluid dynamics to describe energy- 
momentum conservation are written as 

T°£ = -Q„u a , (26) 

where Q„ is the total neutrino emissivity of all processes 
(including the non-equilibrium modified Urea process de- 
scribed by Eq. JDJ); T af) is the energy-momentum tensor 
to be written as (see, e.g., Weinberg 1971): 



rpO 



= Pg 010 + (P + p)u° 



-AT 



a/3 



\rpaf3 



a/3 



AT 



AT j 

cond 



f3S 



-riH a ~<H 



U 7 :S + US : 



9iS 



+ H 01 u c 



) (T ;7 + Tu l; su ). 



(27) 
,(28) 
(29) 



where g a> is the metric tensor, rj is the shear viscosity coeffi- 
cient, k is the thermal conductivity, and H al3 = g al3 + u a vP 
is the projection matrix. In this paper we use r\ — r/ c , where 
the electron shear viscosity r/ c in the stellar core is taken 
from Chugunov & Yakovlev (2005). We neglect the shear vis- 
cosity of neutrons (the proton shear viscosity is even smaller, 
see Flowers & Itoh 1979) because it strongly depends on 
the nuclear interaction model and many-body theory em- 
ployed. A similar problem for heat conduction was discussed 
by Baiko, Haensel & Yakovlev (2001). The neutron shear vis- 
cosity is comparable to the electron shear viscosity (Flowers 
& Itoh 1979), but it cannot change our results qualitatively. 

The use of Eqs. JTSJ, JT7J and (HHJ together with the 
second law of thermodynamics (dp = fj, n dn n + j-i P dn p + 
H e dn a + TdS) can yield the continuity equation for the en- 
tropy in the neutron star core (see, e.g., Landau & Lifshitz 
1959; Weinberg 1971): 



(Su a )., a 

Qbulk 
Qcond 



= (Qbulk + Qshear + Qcond — Qv) /feT, 

= ShAT, Q shcar = (at; 



shear "«> 



( AT cond);/3'"o 



(30) 
(31) 
(32) 



Here 5* is the entropy density and Qbulk is the pulsation 
energy dissipating into heat per unit volume per unit time 
owing to the non-equilibrium modified Urea process. The 
latter term can be interpreted as viscous dissipation due 
to an effective bulk viscosity. We will show below that at 
Sfi <C fcaT it coincides with the term commonly considered 
by other authors (see, e.g., Sawyer 1989 or Haensel et al. 
2001). The term Q s h car describes the dissipation of pulsa- 
tion energy into heat due to the shear viscosity. The term 
Qcond is generally responsible for heat diffusion in the star 
bulk and for the dissipation of pulsation energy due to heat 
conduction. Finally, Q„ is the neutrino emissivity. In this 
work, the quantity Q con d was calculated using an unper- 
turbed metric (the metric Q with v = i/q and A = Ao) 
neglecting temperature variations over a pulsation period. 
The result coincides with the similar expression well-known 
in the cooling theory of non-pulsating neutron stars (see, 
e.g., Thorne 1977; van Riper 1991). These assumptions are 
quite reasonable in the case of a strongly degenerate matter. 
The damping due to heat conduction has been analyzed by 
Cutler & Lindblom (1987) for a more general case of non- 
radial pulsations. The conclusion made by these authors is 
that the contribution of heat conduction to the dissipation 
of pulsation energy can be neglected. 



For the spherically symmetric metric Q, the left-hand 
side of Eq. 1301 can be rewritten as 



d(Se 



A/2x 



M 



1_ d_ 

r 2 dr 



r 2 Se A/2^) 

at 



-(A+i/)/2 



(33) 



In our further treatment, we will use the isothermal ap- 
proximation, in which the redshifted internal temperature 
is taken to be constant over the star bulk: T = Te^ 2 = 
const. This approximation works well for a cooling star of 
the age t > (10 — 50) yrs (see, e.g., Yakovlev et al. 2001; 
Yakovlev & Pethick 2004). The isothermal approximation 
considerably simplifies calculations without any significant 
loss of accuracy (at least for the case of non-equilibrium 
modified Urea processes). Using Eq. 1331 and the integral 
form of Eq. 1301 , and averaging over a pulsation period, we 
arrive at the thermal balance equation: 



d£th _ „ dT 
~dT =CT ~dT 

Ct= I c T dV, 



L phot = AtvR 2 aT s 4 e- o(fl) 



-^phot Li* + W b uik + W shcar , (34) 

(35) 

L v = [~Q„ e"° dV, (36) 



Wbull 



e° dV, W a 



e° dV. 



(37) 



Here, T p h t an d L v are the redshifted photon and neutrino 
luminosities of the star; Wbuik and Wshear denote the heat 
released in the star per unit time owing to the bulk and shear 
viscosities, respectively; a is the Stefan-Boltzmann constant; 
T s is the effective surface temperature (the relationship be- 
tween the surface and internal temperatures is reviewed, 
e.g., by Yakovlev & Pethick 2004). The upper horizontal 
line denotes averaging over a pulsation period. In deriving 
Eq. 13411 . we neglected the terms of the order of ~ T 2 e 2 , 
as compared to those of ~ T 2 , in the left-hand side of the 
equality. The term Q con d in Eq. 1301 leads to the appearance 
of Tphot in Eq. jg). 

The rate of the heat release due to the shear viscosity 



Qshe 



-2r 



d*t(r,t) 



drdt 



2 dt{r,t) r d£{r,t) du 



dt 



6 



u k C k {t) 



dt 
du 



dr 



-2r^- + 2£ fc + r£ fc , 
dr dr 



(38) 



Among the processes contributing to the neutrino emissiv- 
ity Q„, the only process, whose emissivity Q none q can vary 
dramatically over a pulsation period, is the modified Urea 
process (we assume that the direct Urea process is forbid- 
den). The expression for Q noncq is obtained from Eq. i\'2l 
by averaging over the pulsation period P — 2n/u}k with al- 
lowance for y(r,t) — yo cos(u>kt), where yo is a slowly vary- 
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ing function of time: 



110107r 2 yg 
11513 



85057T 4 yt 5257r b yg 315^ y 8 



46052 



46052 



1473664 



(39) 



yo 



C k {t) d5[i(n h0 ,x e0 ) 



"K 2 k-BT 



n-bo ■ 



s^ 2 8 



dr 



r 2 e -,o/2 



(40) 

Similarly, Eqs. (I1H and 1251 result in the expression for the 
heating rate produced by the dissipation of the pulsation 
energy due to deviations from beta-equilibrium: 

_ 14680^ 

Vbulk — 11513 q 

yi 5677T 2 y$ 1057r 4 yo 21tt 6 ^ \ 
2 2936 5872 46976 / ' 1 j 

Analogous expressions for the non-equilibrium direct Urea 
process are presented in the Appendix. 

The quantities Q noncq (yo) and Q bulk (yo) for the mod- 
ified Urea process were calculated numerically by Finzi & 
Wolf (1968). Their results (their Fig. 1) at yo < 1 are correct 
only qualitatively (Haensel 1992) , although they become ex- 
act in the limit of yo S> 1. 

Using Eqs. and EH . one can easily find the neu- 
trino emissivity and the viscous dissipation rate of the non- 
equilibrium modified Urea process for both subthermal or 
suprathermal pulsations (if they are small, i.e. £ <C 1). The 
typical value y of the parameter yo in the stellar core can 
be estimated for the fundamental mode as: 



y ~ 100eC fc (t)/T 9 . 



(42) 



Thus, we have derived the equation describing the ther- 
mal evolution of a pulsating neutron star. This equation de- 
pends on the current pulsation amplitude Ck it) and, hence, 
on the pulsation energy E pu \ B (t) (see Eq. 0). Let us now 
derive the equation to describe the evolution of the pulsa- 
tion energy. In principle, it can be obtained from the "pul- 
sation equation" (59) of Chandrasekhar (1964) by taking 
into account the dissipation terms and considering them as 
small perturbations (generally non-linear in We have 
performed this derivation, but here we will present a much 
simpler derivation following from energy conservation law. 
One should bear in mind that the pulsation energy dis- 
sipates due to the bulk and shear viscosities and is fully 
spent to heat the star. The corresponding terms have al- 
ready been found for the thermal balance equation I34H . 
The same terms, but with the opposite sign, should be valid 
for the damping equation which can be thus presented in 
the form: 

di5p U i s 



At 



-Wbull 



(43) 



The set of Eqs. I|34|l and il l-'il has to be solved to obtain 
self-consistent solutions for the pulsation amplitude Ck and 
temperature T as a function of time. 

Similar equations were formulated, analyzed and solved 
numerically by Finzi & Wolf (1968) under some simplified 
assumptions. In particular, the authors neglected the pul- 
sational damping due to the shear viscosity {Wshe&x = 0). 
They used approximate expressions for L„ and Wbuik (see 
above) and neglected the effects of General Relativity. In 



addition, they used simplified models of neutron stars and 
stellar oscillations. However, their approach was quite suffi- 
cient to understand the main effects of the non-equilibrium 
modified Urea process on the thermal evolution of neutron 
stars and damping of their vibrations. We extend this con- 
sideration using the updated microphysics input with the 
proper treatment of the effects of the shear viscosity and 
General Relativity. 

We can generally write: 

C T = 10 39 a c To erg K" 1 , L„ = 10 40 a L f 9 8 erg s"\ 

£ P ui s = 10 53 a P u\ e 2 Ct erg, E th = C T f/2, 

L v = L v0 (l + ai yl + a,2 yl + 03 y% + «4 yo) > 

W shcar = 10 38 a s u'l yl erg s _1 , 

,„ t f 2 _2 . 4 _ 4 _ 6 8 _ g 

Wbuik = L v0 I - ai y + g 02 y + 2a 3 y + - a 4 y 



y = 10 eCk/Tg. 



(44) 



Here Tg — T/(10 9 K); L v o is the neutrino luminosity of a 
non-pulsating star; ac, aL, ap, as, 01,. . . ,04 are dimension- 
less factors which depend on a stellar model and on a pulsa- 
tion mode. For our model of a neutron star with M = 1.4 Mq 
(the equation of state of Heiselberg & Hjorth- Jensen 1999) 
and for the fundamental pulsation mode we have obtained 
a c = 1.88, a L = 5.34, a P = 1.81, a s = 4.75, oi = 9.18, 
a 2 = 20.0, a 3 = 15.0, a 4 = 3.61. 



5 ANALYTICAL SOLUTIONS AND LIMITING 

CASES 

Before presenting numerical solutions to Eqs. I34H and 14311 . 
let us point out general properties of the solutions and con- 
sider the limiting cases. Numerical values will be given for 
the fundamental mode and for the above model of the star 
with M = 1.4 Mq. 

Equation 1431 describes the damping of pulsations due 
to the bulk and shear viscosities. In the present problem, 
there are no instabilities that could amplify stellar pulsa- 
tions. In contrast, the thermal balance equation 134H permits 
both the stellar cooling (at L„ + Lp hot > Wburk+Wshear) and 
the heating due to the viscous dissipation of the pulsation 
energy (at L v + L pho t < IFbuik + lUshcar). 

One may expect qualitatively different solutions in the 
subthermal (y Q <g 1) and suprathermal (y ^> 1) regimes. 
For the former, we have -E pu i s <C Eth, while for the latter 
-Epuis S> Eth. 



5.1 The modified Urea regime 

A sufficiently hot star has L pho t <§C L u and W B hem <S Wtmik- 
Then the evolution of pulsations and the thermal evolution 
of the star are totally determined by the non-equilibrium 
modified Urea process. The main features of this modi- 
fied Urea regime were analyzed by Finzi & Wolf (1968). 
We present this analysis using more accurate approach (see 
above). 

This regime is conveniently studied by analyzing the 
evolution of T(t) and y (t). Neglecting Lphot and VUshcar, we 
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can rewrite Eqs. 1341 and 1431 as: 



2-E'th dy 

2E th df 
f dt 



L u - W huik 1 + 



-Epuls 

= -L v + W>uik = f s B(y 



T 8 A(y ),(45) 



(46) 



where E tb /E pu \ B = ac/(20 y\ ttp wj); the functions v4(y ) 
and B(y ) are independent of T. Their exact form 
is easily deduced from Eq. (1441 . One immediately has 
dln|/ /dlnT = j4(y )/B(y ), which allows one (in princi- 
ple) to obtain the relation between y Q and T in an integral 
form. 

Equations 1451 and 1461 have two special solutions. 
The first solution is obvious and refers to an ordinary 
non- vibrating (y (t) = 0) cooling neutron star. In this case 



T(t) 



T(0)/ l + 6/%T 9 6 (0)t 



1/6 



(47) 



where /3o = cel /(10 8 ac). We have /3o ~ 1/(1.12 yr), for our 
neutron star model. 

The second solution is realized (Finzi & Wolf 1968) if 
the initial value y Q (0) satisfies the equation 



WbulJ 



1 + 



-Eth 
-Epuls 



1 + 



20 yo o-p w 4 



(48) 



at which A(y ) — and dy(t)/dt — 0. In this case j/ (t) 
remains constant during the modified Urea stage. We will 
denote this specific value of y by y 0L ; it is equal to y 0L « 
0.607 for our neutron star model. In this limiting case T(t) 
and Ck(t) are easily obtained from Eq. 13411 : 



T(t) = T(0)/ (l + 6,ST 9 6 (0)i 
C*(t) = f 9 (t)y 0i /10 2 e, 



1/6 



(49) 
(50) 



where /3 = a h (L u - WWk)/(10 8 a c L v0 ) (» l/(3.05yr), for 
our model). Thus, the internal stellar temperature T(t) and 
the pulsation amplitude Ck (t) simultaneously decrease with 
time, leaving the suprathermality level constant, interme- 
diate between the subthermal and suprathermal pulsation 
regimes. The decrease is power-law (non-exponential). 

The thermal evolution of neutron stars in the two lim- 
iting cases isj-emarkably similar. If the star was born suffi- 
ciently hot (T(0) > 10 9 K) then in a few years after the birth 
the initial temperature becomes forgotten. For the non- 
vibrating star from Eq. we have T 9 (1) (f) w (6/3o t)~ 1/6 , 
while for the vibrating star from Eq. 14911 we have T 9 2 ' (t) w 
(6/3 1) -1 ^ 6 . Thus, the vibrating star stays somewhat hotter, 
?W(t)/fW(t) = (/3o//3) 1/6 (« 1.18 for our model). 

Once the two limiting solutions are obtained, all other 
solutions for the modified Urea regime become clear. If 

I7o(°) > %l> then A (Vo(°)) < and wiU tend to Vol 
from above. If j7 (0) < y 0L , then A(y o (0)) > and y (t) will 
tend to y 0L from below. After y (t) comes sufficiently close 
to y 0L , the stellar evolution is approximately described by 
the limiting solution given by Eqs. I|49|l and l|5()|l . Therefore, 
this limiting solution describes the universal asymptotic be- 
havior of all vibrating neutron stars. 



5.2 The damping of oscillations by the shear 
viscosity 

One may expect qualitatively different solutions for the 
damping due to the bulk viscosity (Wtmik 3> W s hea.t, a 
hot star) and the shear viscosity (VKshoar 2> VKbuik, a cold 
star). Using Eqs. 1441 . it is possible to show that the tem- 
perature T v i sc separating these two regimes (and obeying 
the condition Wbuik ~ Wshear) is approximately equal to 
T v i sc ~7x 10 8 /(1 + y 2 f ls K. For the regime of damping 
due to shear viscosity (T <C Tvisc), Eq. 1431 reduces to a 
linear equation for C'k(t), irrespectively of the value of y : 

dC k (t) 



dt 



C k (t), 



(51) 



where a s hcar ~ 3 x 10 -11 s _1 ~ 1/(1000 yrs) is a constant 
factor. The solution to this equation shows an exponential 
damping, which is independent of y : 



C k (t) = C k (t ) exp 



t^shcar 



dt' 

to fi(t') 



(52) 



5.3 Subthermal Pulsations 

In this case (yo <S 1), Eq. 1411 can be reduced to 



Qbulk — 



14680tt^ 
11513 



= c 



c = 



1 A 

r 2 dr 



-"0/2 



dt 



14680 Q c 



2 

n b o 



dSfi(n b0 , £co) 



r) 



(53) 



(54) 



115137T 2 (k B T) 2 (w fc e-"o/2)2 

The quantity £ can be treated as the bulk viscosity. Equation 
1541 coincides with the corresponding expression of Sawyer 
(1989) and Haensel et al. (2001). If the temperature remains 
constant during the damping, Eqs. J7J and 1431 yield an 
exponential fall of the pulsation amplitude Ck(t), which is 
often discussed in literature (see, e.g., Cutler et al. 1990). 

According to Eq. 1441 . the subthermal regime is char- 
acterized by L v w L„o 3> Wbulk- In this case pulsations do 
not affect the neutrino luminosity, and the energy dissipation 
due to the bulk viscosity cannot produce a considerable stel- 
lar heating. The dissipation due to the shear viscosity is also 
too weak, IU s hcar L v q. For these reasons, the pulsations 
do not change significantly the thermal balance equation 
1341 and the thermal evolution of the star. At the neutrino 
cooling stage (when L„o 3> iphot, which happens at t < 10 
yrs) we get the well-known formula 1471 for non-superfiuid 
neutron stars that cool via the modified Urea process. It can 
be rewritten as (see, e.g., Yakovlev & Pethick 2004) 



t = CtT/(6L„ ) ~ 1 yr/T 9 6 



(55) 



This value of t can be considered as a characteristic cooling 
time tcooi of the star with the internal temperature T. For a 
hot star with Wbuik 3> WWar (T 2> T v i sc ; the modified Urea 
regime), Eq. 1431 gives the characteristic pulsation damping 
time 



ipuls ~ -Epuls/Wbulk ~ tcool- 



(56) 



Therefore, the internal temperature T and the typical dis- 
balance of the chemical potentials 8fi decrease with approxi- 
mately the same characteristic time t coo i (see, e.g., Yakovlev 
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et al. 2001). The parameter y oc 5fi/T, which describes the 
"level" of pulsations relative to the thermal "level", should 
tend to the limiting value y 0L (Sect. I5TT1 . The damping of 
pulsations obeys the power law (rather than exponential, as 
would be in the absence of the thermal evolution). This is 
because the viscous damping rate strongly depends on tem- 
perature, Wbuik oc T 6 ^ 

In a cooler star (T < T v isc, Wshcar > Wbuik), the damp- 
ing of subthermal pulsations is due to the shear viscosity 
and occurs, according to Eq. 1521 . more abruptly (exponen- 
tially), decreasing the pulsation level y . 

5.4 Suprathermal Pulsations 

In this case, the quantity Q bulk cannot be generally de- 
scribed by an expression of the type of Eq. I|53|l . Strictly 
speaking, we cannot introduce a bulk viscosity f, but Eq. 
14111 adequately describes the rate of the pulsation energy 
dissipation due to the modified Urea process. Nevertheless, 
at least for radial suprathermal pulsations, the quantity 
Qbuik can t> e formally calculated from Eq. 1531 . as before, 
with the effective bulk viscosity £ given by Eq. 1541 with an 
additional factor Qbuik/Qbuik(j7o ~~ * 0)- ^ n tne suprathermal 
regime, the effective bulk viscosity and the viscous dissi- 
pation rate appear to be much larger than in the subther- 
mal regime, as was pointed out by Haensel, Levenfish & 
Yakovlev (2002) . However, there is an omission in their Eqs. 
(13)-(15) for the effective bulk viscosity in the suprathermal 
regime: the authors should have introduced an additional 
factor ~ (1 + ?7q). This does not change qualitatively their 
principal results. Nevertheless, we stress that while analyz- 
ing the damping of pulsations, one should account for the 
thermal evolution of the star. Accordingly, in the suprather- 
mal regime Eq. (16) of Haensel et al. (2002) gives the char- 
acteristic time of non-exponential damping. 

The pulsation equation 143H in the suprathermal regime 
(£ puls > E th , y Q > 1) at W bu lk > W s hear (f > T visc ; the 
modified Urea regime) can be rewritten as 

— ~tt — = -Q(buikC fc (t), (57) 

where Qbuik ~ 3 x 10 4 £ 6 /u>| s _1 is a constant factor. As- 
suming Cfc(0) = 1 we get: 

C fc (t) « (l + 3Q bulk t)" 1/6 . (58) 

This solution describes a slow (power law) fall of the 
pulsation amplitude ~ t -1 / 6 with the characteristic time 
l/(3Qbuik). In this regime, the stellar heating always dom- 
inates over the cooling, with the heating rate Wbuik ~ 
8/3 L v oc 5/j, nearly independent of temperature (being de- 
termined by the disbalance of the chemical potentials Sfj,). 
This result was obtained by Finzi & Wolf (1968). The power 
law decrease of Ck{i) is associated with a strong depen- 
dence of Wbuik on S/j, (which mimics the dependence on T 
in the subthermal regime). The relative pulsation amplitude 
should decrease, and the star should evolve to the subther- 
mal regime (y — > y 0L ; Sect. 15. lV 

In a rather cool star, the damping of pulsations due to 
the shear viscosity dominates over the damping caused by 
the bulk viscosity (Wshear 3> Wbuik, T < T v i sc ). The shear 
viscous damping is exponential, according to Eq. 1511 . so 
that the star rapidly evolves to the subthermal regime. 



6 RESULTS 

Generally, the set of Eqs. I34H and 1431 has no analytical 
solution, and we have to solve it numerically. We have mod- 
ified the isothermal version of our cooling code (for details, 
see the review of Yakovlev et al. 1999) by including a block 
for solving the damping equation (1438 . Our code calculates 
the stellar surface temperature T a °° (redshifted for a distant 
observer), as a function of time t, as well as Ck(t). All the 
computations presented in this section are for the funda- 
mental mode of radial pulsations. Computations for higher 
modes will not lead to qualitatively different conclusions. 

The left panel of Fig. 2 shows the thermal evolution 
paths of a neutron star (M = 1.4 Mg), which differ in the 
initial internal temperature T(0) = To and the initial rela- 
tive amplitude of pulsations e (see Eq. J3J). The right panel 
presents Cfe(i) curves for the same models. The dotted curve 
on the left panel shows the cooling of a non-pulsating star 
(in the isothermal approximation). The circle indicates the 
observations of the Vela pulsar. References to the original 
observations can be found in Gusakov et al. (2004). 

The solid lines in both panels of Fig. 2 are for the model 
with Tb = 10 10 K and e = 0.01. This model describes a star 
which was born hot and strongly pulsating. The initial pul- 
sation energy is about twice lower than its initial thermal 
energy, and the star is in an intermediate pulsation regime, 
between the supra- and subthermal regimes, with S/j, ~ &bT. 
The heating due to viscous dissipation is not as fast as the 
neutrino cooling due to the non-equilibrium modified Urea 
process, and the star is cooling down. The main contribu- 
tion to the dissipation at the initial stage is produced by the 
bulk viscosity. The maximum difference between the surface 
temperatures of such star and a non-pulsating star occurs 
at t < 1000 yrs. During this period of time, Sfi remains of 
the order of k^T (y w Vol)- At t J> 1000 yrs, the damp- 
ing begins to be determined by the shear viscosity, which 
is not so temperature-dependent as the bulk viscosity. This 
leads to the exponential damping in the subthermal regime 
(i? P uis/-Bth < 1; see the right panel of Fig. 2). 

The short dashed lines in Fig. 2 correspond to an ini- 
tially cold star with To = 10 6 K and e = 0.01. The ini- 
tial ratio of the pulsation energy to the thermal energy is 
Spuiso/Stho ~ 5 x 10 7 , i.e., the star pulsates in a strongly 
suprathermal regime. As follows from Eq. 144L at low tem- 
peratures we have Wbuik "C Wshear, and the star is initially 
heated up by the shear viscosity. The heating due to the 
bulk viscosity starts to dominate only at T > 5 x 10 s K. 
After heating to T w 1.7 x 10 9 K in t ~ 1 month, the star 
appears in the intermediate regime with Sfi ~ feT and be- 
gins to cool down. At t > 10 yrs, the star starts evolving 
along the same "universal" path (y m y 0L ) as in the first 
model. 

The long dashed curves are obtained for the same ini- 
tial conditions but in the "naive" approximation neglecting 
non-linear effects in non-equilibrium beta-processes. In par- 
ticular, the neutrino luminosity is taken to be L v = L„o 
and the damping due to the bulk viscosity is determined 
by Eqs. 153H and I54H . One can see that this approximation 
leads to qualitatively incorrect results. The viscous heating 
during the first year after the pulsation excitation is much 
slower than in the scenario with non-linear effects, and the 
star heats up slowly. The neutrino luminosity is also lower, 
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Figure 2. Temporal evolution of the surface temperature (left) and pulsation amplitude (right) for three initial conditions: To = 10 10 
K, e = 0.01 (solid lines); To = 10 6 K, e = 0.01 (dashed lines); To = 10 6 K, e = 0.0001 (dash-dotted lines). The short dashed lines denote 
a self-consistent calculation, the long dashed lines are calculated neglecting non-linear effects due to deviations from the beta-equilibrium 
(see text). The dotted line on the left panel shows the cooling of a non-pulsating star. The full circle denotes the observations of the Vela 
pulsar. The insert displays the pulsation damping at t < 30 yrs in more details. 



which enables the star to heat to higher temperatures. In 
fact, the slow damping due to the bulk viscosity does not 
decrease the pulsation amplitude during the first year. 

The dash-and-dot lines in Fig. 2 refer to the cold star 
with f = 10 6 K and e = 0.0001. The initial ratio of the pul- 
sation and thermal energies is -Epuiso/^tho ~ 5 x 10 3 , which 
means that the star is initially in the suprathermal regime. 
Nevertheless, the pulsation energy -B pu i s o ~ 5 x 10 45 erg is 
insufficient to heat the star to a temperature at which the 
damping is determined by the bulk viscosity. For this rea- 
son, the pulsation energy is damped by the shear viscosity. 
The damping of pulsations takes ~100 years (see the right 
panel of Fig. 2). At t > 100 yrs the star cools via photon 
emission from the surface. It is clear from the left panel that 
this model can, in principle, explain the surface temperature 
of a neutron star with the same thermal X-ray luminosity 
as the Vela pulsar but with different history. For example, 
it may be an old and cold isolated neutron star, in which 
radial pulsations have been excited. In ~10 years after the 
excitation, the star will be heated to the temperature of the 
Vela pulsar. In 100 years, the pulsations will die out but the 
star will stay warm for ~ 10 years before it starts cooling 
down noticeably. It should be emphasized that these results 
will not change if we take a lower initial temperature, e.g., 
f = 10 4 K. This star will also acquire the surface tem- 
perature T s °° ~ 7 x 10 5 K in a year and will emit in soft 
X-rays. 

For a correct calculation of the pulsation damping, one 
should take into account the thermal evolution of the star. 
The evolutionary effects are especially important when the 
damping is determined by non-equilibrium beta-processes. 



They are relatively weak only in the subthermal regime, pro- 
vided the damping is produced by shear viscosity. 

These statements are also illustrated in Fig. 3 which 
shows the pulsation damping for a star with the initial in- 
ternal temperature T(0) = 10 9 K and the initial relative pul- 
sation amplitude e = 0.01. The initial ratio of the pulsation- 
to-thermal energy is -B pu i s o/-Btho ~ 50, indicating that the 
star is pulsating in a slightly suprathermal regime. The solid 
line is the result of a self-consistent solution of the thermal 
evolution and damping equations. The damping is power 
law for about 100 years; afterwards the damping is deter- 
mined by the shear viscosity and becomes exponential. The 
pulsations die out completely in ~1000 years. 

The dotted line in Fig. 3 shows the solution to the 
damping equation neglecting the thermal evolution, at a 
constant internal temperature T = T(0). In this case, the 
pulsations are first damped by the bulk viscosity and then 
by the shear viscosity in ~ 30 years. 

The short dashed curve is obtained by taking into ac- 
count the thermal evolution and damping, but neglecting 
the non-linear effects in non-equilibrium beta-processes. For 
about 100 years, the damping is governed by the bulk vis- 
cosity; it is power law, but slower than with the non-linear 
effects. Later, the shear viscosity becomes important, lead- 
ing to the exponential damping, nearly the same as with the 
non- linear effects. 

Finally, the long dashed curve is calculated neglecting 
both the thermal evolution and the non-linear effects. Like in 
the case with these effects (the dotted curve), the damping 
is steep (exponential), taking about 30 years, but occurs 
slightly slower (the long dashed curve is above the dotted 
curve). 
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Figure 3. Temporal evolution of the pulsation amplitude for the initial conditions: Tq = 10 K, e = 0.01. Solid line presents a consistent 
calculation with account for the thermal evolution and non-linear deviations from beta-equilibrium (see text); dotted line shows the same 
but the internal stellar temperature is fixed; short dashed line is the same as the solid line but without non-linear effects; long dashed 
line is obtained neglecting non-linear effects at the fixed temperature. The insert displays the pulsation damping at t <J 30 yrs in more 
details. 



7 SUMMARY 

Extending the consideration of Finzi & Wolf (1968) we 
have analyzed the thermal evolution of a non-superfluid star 
which undergoes small-amplitude radial pulsations. We have 
derived a set of equations to describe the thermal evolution 
and the damping of pulsations in the frame of General Rela- 
tivity. We have included the effects of non-linear deviations 
from beta-equilibrium in the modified Urea process on the 
neutrino luminosity and on the pulsation energy dissipation 
due to the bulk viscosity in the stellar core. We have also 
taken into account the dissipation due to the shear viscos- 
ity and the associated heating. A set of equations for the 
evolution of a neutron star with a nucleon core, in which 
the direct Urea process is forbidden, has been analyzed and 
solved analytically and numerically. 

We have shown that the evolution of a pulsating star 
strongly depends on the degree of non-linearity of the non- 
equilibrium modified Urea process and on the nature of the 
pulsation damping (the shear or bulk viscosity) . In the non- 
linear regime, the star may be considerably heated by the 
pulsation energy dissipation but it is always cooled down in 
the linear regime. Pulsations of a hot star are damped by the 
bulk viscosity in both the linear and non-linear regimes, and 



this process is rather slow (power law). In a cooler star, the 
damping is produced by the shear viscosity and goes much 
faster (exponentially). The characteristic times of damping 
of the fundamental mode lie within 100-1000 years. 

We have not discussed here the specific damping mech- 
anism via the ambipolar diffusion of electrons and protons 
relative to neutrons when the averaged (over the period) 
chemical composition of the stellar matter remains constant 
in time. As far as we know, this mechanism of the pulsation 
damping has not been analyzed in the literature. However, 
it may be as efficient as the damping by the shear viscosity, 
at least in the suprathermal regime. We will consider this 
problem in a separate publication. 

The analysis presented here is based on a simplified 
model. In particular, if the direct Urea process is open in the 
stellar core or if the core contains hyperons or quarks, the 
bulk viscosity can be many orders of magnitude higher than 
discussed here (see, e.g., Haensel et al. 2002 and references 
therein). The results may also differ significantly for super- 
fluid neutron stars, because superfluidity drastically changes 
the reaction rates in dense matter and, hence, its kinetic 
properties, including the viscosity. It would also be instruc- 
tive to consider other types of neutron star pulsations, pri- 
marily r-modes. They can be accompanied by the emission of 
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gravitational waves (see, e.g., Andersson & Kokkotas 2001) 
which can, in principle, be registered by gravitational detec- 
tors of new generation. We expect to continue the analysis 
of the evolution of pulsating neutron stars. 
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APPENDIX A: NON-EQUILIBRIUM DIRECT 
URCA PROCESS 

If the non-equilibrium direct Urea process is allowed in a pul- 
sating neutron star, it can also be described by the quantities 
Qnoneq and Qbuik given by Eqs. H12H and H31L respectively. 
These quantities are taken from Yakovlev et al. (2001) and 
are denoted here as Qnoneq(jz) and Q b ui k (y)- The averaging 
over a pulsation period yields 
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where Qcq' is the neutrino emissivity of the direct Urea 
process and (as before) yo = 8 /i / (n 2 k^T) . In a pulsating 

star with the allowed direct Urea process, and Q^ulk 

should be included into the quantities Q lloncq and Q bulk 
given by Eqs. 11341 and 1431 . respectively. In the absence 
of nucleon superfluidity, the contribution of the direct Urea 
process into Q^ ncq and Q^uik is 5 - 7 orders of magnitude 
greater than that of the modified Urea process. 
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